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We study a vectorial gauge theory with gauge group SU(iV c ) and a variable number, Nf, of 
massless fermions in the fundamental representation of this group. Using approximate solutions 
of Schwinger-Dyson and Bethe-Salpeter equations, we calculate meson masses and investigate how 
these depend on Nf. We focus on the range of Nf extending from near the boundary with a non- 
Abelian Coulomb phase, where the theory exhibits a slowly running ("walking") gauge coupling, 
toward smaller values where the theory has non-walking behavior. Our results include determi- 
nations of the masses of the lowest-lying flavor-adjoint mesons with J PC = CT+, I -- , ++ , and 
1 ++ (the generalized it, p, ao, and ai). Related results are given for flavor-singlet mesons and for 
the generalization of f^. These results give insight into the change from walking to non- walking 
behavior in a gauge theory, as a function of Nf. 
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I. INTRODUCTION 

We consider a (3 + l)-dimensional vectorial gauge the- 
ory (at zero temperature and chemical potential) with 
the gauge group SU(A^ C ) and Nf massless fermions trans- 
forming according to the fundamental representation of 
this group. For N c — 3, if one took Nf = 2, this 
would be an approximation to actual quantum chromo- 
dynamics (QCD) with just the u and d quarks, since 
their current quark masses are small compared with the 
scale Aqcd — 400 MeV. We restrict here to the range 
Nt < (U/2)N C for which the theory is asymptotically 
free. An analysis using the two-loop beta function and 
Schwinger-Dyson equation (reviewed below) leads to the 
inference that for Nf in this range, the theory includes 
two phases: (i) for < Nf < Nf >cr a phase with confine- 
ment and spontaneous chiral symmetry breaking (S^SB); 
and (ii) for Nf jCr < N f < (U/2)N C a non-Abelian 
Coulomb phase with no confinement or spontaneous chi- 
ral symmetry breaking. We shall refer to Nf^ cr , the crit- 
ical value of Nf, as the boundary of the non-Abelian 
Coulomb (conformal) phase Q . Here we take electroweak 
interactions to be turned off. We denote the fermions as 
ff with a = 1,...,N C and i = 1, ...,Nf. The theory has 
an SXJ(N f ) L x SV(Nf) R x U(l)y global symmetry (the 
U(1)a being explicitly broken by instantons), which is 
spontaneously broken to SU(7V/)y x U(l)y by the for- 
mation of a bilinear fermion condensate. 

For Nf slightly less than Nf >cr , the theory exhibits 
an approximate infrared (IR) fixed point with resultant 
walking behavior. That is, as the energy scale /i decreases 
from large values, a = g 2 /(4ir) (g being the SU(iV c ) 
gauge coupling) grows to be 0(1) at a scale A, but in- 
creases only rather slowly as fj, decreases below this scale, 
so that there is an extended interval in energy below A 
where a is large, but slowly varying. Associated with 
this slowly running behavior, the resultant dynamically 



generated fermion mass, £, is much smaller than A. In 
addition to its intrinsic field-theoretic interest, this walk- 
ing behavior has played an important role in theories of 
dynamical electroweak symmetry breaking 0-0- As Nf 
approaches Nf tCr from below, quantities with dimensions 
of mass vanish continuously; i.e., the chiral phase transi- 
tion separating phases (i) and (ii) is continous. Recently, 
meson masses and other quantities such as the general- 
ized pseudoscalar decay constant f n were calculated in 
the walking limit of an SU(iV c ) gauge theory §. 

In the present paper we shall investigate how meson 
masses and other quantities change as one decreases Nf 
below Nf jCr , moving away from the boundary, as a func- 
tion of Nf, between phases (i) and (ii), deeper into the 
confined phase. Our paper is thus a study of the change 
(crossover) between the walking behavior that occurs 
near to this boundary, and the non-walking behavior that 
occurs for smaller Nf. In a non- walking (asymptotically 
free, confining) theory such as real QCD, as the energy 
scale \i decreases through A, a increases rapidly through 
values of order unity, triggering spontaneous chiral sym- 
metry breaking on this scale, so that £ ~ A. This is quite 
different from a theory with walking, in which £ « A. 
Our basic calculational methods are essentially the same 
as those employed in Ref. |8j, i.e., we use the Schwinger- 
Dyson (SD) equation to compute the dynamical fermion 
mass £ (generalized constituent quark mass) and then 
insert this into the Bethe-Salpeter (BS) equation to ob- 
tain the masses of the low-lying mesons. We restrict to 
an interval of Nf values for which the theory has an in- 
frared fixed point (as calculated from the beta function, 
to be discussed further below). The reason for this is 
that it makes our calculations more robust since for our 
interval of Nf we can avoid having to introduce a cutoff 
on the growth of a in the infrared. If one uses Schwinger- 
Dyson and Bethe-Salpeter equations to explore a region 
of Nf where the beta function does not have an infrared 
fixed point, one must use such an IR cutoff, which leads 



to cutoff-dependence of the results. For definiteness, we 
shall take N c — 3; however, as will be seen, N c only enters 
indirectly, via the dependence of the value of the infrared 
fixed point a* (eq. I|2.6[) below) on N c . Hence, our find- 
ings may also be applied in a straightforward way, with 
appropriate changes in the value of a*, to an SU(AT C ) 
gauge theory with a different value of N c . 

We mention some background and related work. Many 
studies have investigated the hadron mass spectrum for 
QCD with Nf = 2 light quarks. Among the earliest were 
static quark models @, and bag models Lattice 
gau ge theory has provided an especially powerful method 
|l2|. The Schwinger-Dyson and Bethe-Salpeter equations 
have been used for many years to study spontaneous chi- 
ral symmetry breaking and relativistic bound states in 
field theories (a partial list of papers and reviews includes 
0- 01; E3"E3)- I n particular, the Bethe-Salpeter equa- 
tion has been used to calculate meson masses in QCD 
HJ-IHl. For the walking limit, in addition to Refs. ||, 
these methods have also been used in connection with 
spectral function sum rules to study the ir + — ir° mass 
difference [3(| and the S parameter (equivalently, the chi- 
ral Lagrangian coefficient Lio) [3l| . 

This paper is organized as follows. In Section II we 
review some background material concerning the beta 
function, approximate infrared fixed point, and walk- 
ing behavior. Section III includes a discussion of the 
Schwinger-Dyson equation and our solution of it, as well 
as our calculation of the pseudoscalar decay constant fp, 
the generalization of /„.. In Section IV we present our cal- 
culation of meson masses using the Bethe-Salpeter equa- 
tion. Section V contains some further remarks and our 
conclusions. 



II. PRELIMINARIES 

In order to study meson masses and other quantities as 
one moves away from the boundary between the confined 
phase with spontaneous chiral symmetry breaking and 
the non-Abelian Coulomb phase, it is first necessary to 
know as accurately as possible where this boundary lies, 
as a function of Nf, i.e., to know the value of Nf >cr . We 
first review the estimate M| based on using the two-loop 
SU(iV c ) beta function [MEI 
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bi 



(2.1) 



where t = In /i, with /i the energy scale. The two terms 
listed are scheme-independent. (Two higher-order terms 
in /? have been calculated but are scheme-dependent; in- 
clusion of these does not significantly affect our results.) 
For the relevant case of an asymptotically free theory, 
bo > so that an infrared fixed point exists if 6i < 0. 
This coefficient b\ is positive for < Nf < NfjR, where 
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(2.2) 



and negative for larger Nf. For iV c = 3, Nfjp ~ 8.1 
[23. The value of a at this IR fixed point, denoted a*, is 
given by a* = —Anbo/bx. Substituting the known values 
of these terms, one has 



-4tt(11JV c - 2N f ) 



347V 2 - l3N c N f + 3Nc X N f 



Solving eq. (|2.3|l for Nf in terms of a* yields 



N f = 



2N? [17N c (a+/n) + 22] 
(13A c 2 -3)(a*/7r) + 8AT c 



(2.3) 



(2.4) 



In the one-gluon exchange approximation, the Schwinger- 
Dyson gap equation for the inverse propagator of a 
fcrmion transforming according to the representation R 
of SU(AT C ) has a nonzero solution for the dynamically 
generated fermion mass, which is an order parameter for 
spontaneous chiral symmetry breaking, if a > a cr , where 
a cr is given by 



3a cr C 2 (R) 



= 1, 



(2.5) 



and C2{R) denotes the quadratic Casimir invariant for 
the representation R j3f|. Using C2(fund.) = C%f — 
(N^ — 1)/ (2N C ) for the fundamental representation yields 



2nN r 



3(7V 2 - 1) ■ 



(2.6) 



For the case N c = 3 that we use for definiteness here, 
eq. (|2.6|l gives a cr — tt/A ~ 0.79. To estimate Nf cr , one 
solves the equation a* = a cr , yielding the result Q 



N S,cr = 



2iV c (50iV c 2 - 33) 
5(5A^ C 2 - 3) 



(2.7) 



For N c — 3 this gives Nr iCr ~ 11.9. These estimates 
are only rough, in view of the strongly coupled nature 
of the physics. Effects of higher-order gluon exchanges 
have been studied in Ref. [jjj]- These calculations are 
scmi-pcrturbative and do not include instanton effects. 
It is known that instantons enhance the formation of 
the bilinear fermion condensates 36] , which suggests that 
their inclusion would expand the phase with confinement 
and spontaneous chiral symmetry breaking, i.e., would 
increase Nf cr somewhat relative to the value obtained 
from the two-loop beta function and gap equation. In 
principle, lattice gauge simulations provide a way to de- 
termine Nf >cr , but the groups that have studied this have 
not reached a consensus [37l|-[39^. 

In our analysis, what we actually vary is the value of 
the approximate IR fixed point a*, which depends para- 
metrically on Nf. Thus, although our SD and BS equa- 
tions are semi-perturbative, the analysis is self-consistent 
in the sense that our a cr really is the value at which, 
in our approximation, one passes from the confinement 
phase to the non-Abelian Coulomb phase, and our values 
of a do span the interval over which there is a crossover 
from walking to QCD-like (i.e., non- walking) behavior. 
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We elaborate here on the origin of the walking behav- 
ior. Since the theory is asymptotically freee, it follows 
that as the energy scale p decreases from values 3> A, 
a increases. If Nf < Nfjn, there is no perturbative 
IR fixed point. If NfjR < Nf < Nf tCr , as the energy 
scale decreases toward zero, the coupling a approaches 
a*, which is larger than a cr - The coupling a* is only an 
approximate IR fixed point since, as a increases past a cr 
and the fermion condensate forms, the fermions gain a 
dynamical mass E so that in the low-energy effective the- 
ory applicable for energy scales fj, < E, one integrates out 
these fermions and is left with a pure gluonic SU(iV c ) the- 
ory with a different beta function, such that a increases 
further. If TV/ > Nf tCr (and smaller than (ll/2)iV c ), the 
theory is in the non-Abelian Coulomb phase and a* is an 
exact IR fixed point. In the case that Nt is only slightly 
less than Nf jCr , or equivalently, a* is only slightly greater 
than a cr , it follows that as the energy scale decreases and 
a approaches a* from below, the rate of increase of a, 
i.e., |/3 1, decreases, so that the theory has a large coupling 
a ~ 0(1) which, however, runs very slowly. 

As is evident from the above results, decreasing Nf 
below Nf t cr has the effect of increasing a* and thus mov- 
ing the theory deeper in the phase with confinement and 
spontaneous chiral symmetry breaking, away from the 
boundary with the non-Abelian Coulomb phase. This is 
the key parametric dependence that we shall use for our 
study. Our aim is to investigate how meson masses and 
other observable quantities depend on Nf in the crossover 
region; operationally, what we actually vary is a*. In 
Ref. [2| the range of a* used for the calculation of me- 
son masses was chosen to be 0.89 < a* < 1.0, an interval 
where there is pronounced walking behavior. For the case 
N c = 3 considered in Ref. and here, given the above- 
mentioned value, a cr = 7r/4, it follows that this lower 
limit, a* = 0.89, is about 12 % greater than this critical 
coupling. The reason for this choice of lower limit on a* 
was that the hadron masses become exponentially small 
relative to the scale A as a* — a cr — ► + , rendering nu- 
merical evaluations of the relevant integrals increasingly 
difficult in this limit. For our study of the shift away 
from walking behavior we consider an interval extending 
to larger couplings, from a* = 1.0 to a* = 2.5. Our 
upper limit is chosen in order for the ladder approxima- 
tion used in our solutions of the Schwinger-Dyson and 
Bethe-Salpeter equations to have reasonable reliability. 
From eq. 1)2.4(1 it follows that a* = 0.89 corresponds to 
Nf = 11.65, about 2 % less than Nf )Cr . For a coupling 
as large as a* = 2.5, the semi-perturbative methods used 
to derive eqs. 1)2.3)1 and (|2.4(l are subject to large cor- 
rections from higher-order perturbative, and from non- 
perturbative, contributions; recognizing this, the above 
upper limit of a* corresponds formally to Nt ~ 9.8, a 
roughly 20 % reduction from Nf :Cr = 11.9. 

Since the chiral transition which occurs as Nf in- 
creases through Nf :Cr is second-order (continuous), and 
since there is no spontaneous chiral symmetry breaking 
in the non-Abelian Coulomb phase, it follows that as 



Nf /* Nf yCr , (i) the masses of all hadron states vanish 
continuously; and (ii) hadron states that are related to 
each other by a parity reflection become degenerate. 



III. SCHWINGER-DYSON EQUATION 

We first use the Schwinger-Dyson equation for the 
fermion propagator to calculate the dynamically gener- 
ated mass S of this fermion. This extends the calcula- 
tion of these quantities in Ref. [8j to smaller Nf and, 
accordingly, larger a*. The inverse fermion propagator 
is Sf(p)^ 1 — A(p 2 )^> — B(p 2 ). We approximate the full 
Schwinger-Dyson equation by using an effective running 
coupling and the lowest-order gluon propagator: 



Sf( P y 
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(2tt)< 



g 2 (p,q)D^(p-q)rSf(q)Y 
(3.1) 

We use the Landau gauge for the gluon propagator 
D^(k), i.e., D^ v = (-g^ + fc M fc„/fc 2 )/fc 2 because this 
simplifies the calculation. The physical results are, of 
course, gauge-invariant (e.g., [Tflj). Equation ()3.1)> yields 
two separate equations for A(p 2 ) and B(p 2 ). As in Ref. 

we make the ansatz for the running coupling, after 
Euclidean rotation, 



a(pE,qE) 



a{p E 



(3.2) 



where the subscript denotes Euclidean. Since a would 
naturally depend on the gluon momentum squared, (p - 



q) 2 



p 2 + q 2 — 2p-q, the functional form 1)3. 2|) amounts to 



dropping the scalar product term, — 2p ■ q. This is a par- 
ticularly reasonable approximation in the case of a walk- 
ing gauge theory because most of the contribution to the 
integral (|3.1|l comes from a region of Euclidean momenta 
where a is nearly constant. Hence, the shift upward or 
downward due to the —2p-q term in the argument of a 
has very little effect on the value of this coupling for the 
range of momenta that make the most important contri- 
bution to the integral. The approximation l|3.2|l enables 
one to carry out the angular integration, obtaining the re- 
sults A(p 2 E ) — 1 and, for B(p 2 E ) = E(p|,), setting x =p% 
and y \ 
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a(x + y) Y,{y) 



(3.3) 



In terms of the momentum-scale-dependent fermion mass 
T,(p E ), we define the dynamical mass E as 



E = E(p| = E 2 ) 



(3.4) 



As noted above, in the walking region, over most of 
the range of integration over qs in eq. I|3.3|l below A, the 
running coupling a is approximately constant and equal 
to its fixed-point value, a* (see Fig. 2 of Ref. @). This 
means that in the walking region one does not have to 
introduce any infrared cutoff on the growth of a, as was 
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necessary in earlier studies of the Schwinger-Dyson and 
Bethe-Salpeter equations for regular QCD [18j , |2l| , [24| • 
In the interval ^ S C A, the fermions decouple, hav- 
ing gained dynamical masses E, and in this low-energy 
theory with the fermions integrated out, the resultant 
a would evolve away from a* as calculated via the per- 
turbative beta function. However, since this dynamical 
mass scale is much smaller than A in a walking theory, 
it follows that this lowest range of the integration over 
qE makes a negligibly small contribution to the entire 
integral. One can thus employ the approximation of us- 
ing the same functional form for a down to qE = in 
the integral. This convenient feature does not hold if Nf 
decreases very far below Nf^ cr , i.e., a* increases too far 
above a cr . 

Having made the approximation of using the same 
functional form for a for ks in the range < < E 
as in the range fc^ > E, and solving for a(fc^) from the 
two-loop beta function, one finds that, in terms of the 
variable ln(fc^/A), it increases rather quickly from small 
values to values of 0(1) as ks decreases below A. This 
motivates an additional simplification, namely approxi- 
mating a{k,E) as the step function, 



t{k E ) = a*6»(A - k E ) 



(3.5) 



Then as a* \ a cr , if one also approximates the denom- 
inator of the fermion propagator in eq. (|3.3|) as (q E + T, 2 ), 
i.e., one sets Yi(q E ) = E in this denominator, then the so- 
lution is 0-0 



E = cA exp 



-1/2 



(3.6) 



where c is a constant. 

One next discretizes the Schwinger-Dyson equation 
and solves it using iterative numerical methods, as de- 
scribed in Ref. p|. In Fig. we show the solution 
for the dynamical fermion mass E as a function of a*. 
A fit to the numerical solution in the walking region 
0.89 < a* < 1.0 in Ref. Q found agreement with the 
functional form l|3.6|) with c = 4.0. Our calculations for 
larger a* show the expected shift away from walking be- 
havior. This shift is evident in Fig. ^for a* larger than 
about 1.2. Note that our solution of the full Schwinger- 
Dyson equation does not make the approximation of set- 
ting Yi{q 2 E ) — E in the fermion propagator denominator 
but instead incorporates the full functional dependence 
of T,(q E ). In real QCD, precision fits to deep inelastic 
lepton scattering data, hadronic decays of the Z, etc. 
probe the theory in momentum regions where Nt = 4 
or Nf — 5, and yield, for the effective ^/-dependent 

280 MeV, with 



scale A.q CD 



200 MeV and A^ D 



larger values for Aq^, with Nf — 3, 2. In actual QCD 

one thus has E/A^' ~ 0(1) for these low values of 
Nf. These contrast with the limiting walking behavior, 
in which E <C A, as indicated in eq. I|3.6[) . Our calcu- 
lation of E, shown in Fig. ^ shows that E/A increases 
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FIG. 1: Numerical solutions for E, for several values of q, 
(indicated by (>)• By way of comparison, we show, as the 
dotted curve, the solution 13.61 with c = 4.0 derived from a 
fit to the results in the interval 0.89 < a* < 1.0. See text for 
further discussion. 



substantially, by about a factor of 30, from a value of 
about 0.01 at a* = 1.0 to 0.32 at a* = 2.5, much closer 
to the value of 0(1) for this ratio in QCD. 

Another quantity of interest is the pseudoscalar de- 
cay constant fp, the A^-flavor generalization of the pion 
decay constant. For Nt = 2 QCD this is defined as 
(0|Jjj|7r fc (g)) = if^S jk where 1 < j, k < 3 are SU(2) 
isospin indices. Here, we use a generalization of this 
definition, with the symbol f„ replaced by fp and the 
SU(Nf) isospin indices in the range 1 < j, k < Nf — 1. In 
QCD, one rough measure of the dynamical (constituent) 
quark mass is E ~ Mn/N c ~ 313 MeV, where M/v is the 



nucleon mass. An alternate definition sets E ~ M p /2\ 
this would yield a somewhat larger value. Here we use 
the estimate E ~ 330 MeV. Using the measured value 
fa ~ 92.4 ± 0.3 MeV [H, one thus has 
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(3.7) 



QCD 



An approximate relation connecting E and fp is |44j 
(with y=k 2 E ) 



5? (y) - Jf £ [£»(»)] 



[y + E2(y) 



(3.8) 



The integration is rendered finite by the softness of the 
dynamical mass E(fc|»), which behaves for fc > E as 



E(*4)ocE( — 
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(3.9) 



where 7 is the anomalous dimension of the bilinear op- 
erator //, having the value 7 ~ 1 in the walking regime 
and decreasing toward zero at very large energy scales 
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FIG. 2: Values of fp calculated from eq. 1)3.8^ for several val- 
ues of a, (indicated by 0). For comparison, we show, as the 
dotted line, the analytic solution given by eq. 113.611 with E re- 
placed by fp and c = 1.5 derived from a fit to the calculations 
for 0.89 < a, < 2.5. See text for further discussion. 



fj, » A (since 7 is a power series in a and a — > in 
this limit due to the asymptotic freedom of the theory) . 
Thus, the relation l|3.8[l suggests that for QCD 
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For N c = 3, this is £//„ — 2vt/v / 3 ~ 3.6, which agrees, 
to within the theoretical uncertainties, with experiment. 
In QCD, with A( 2 ) ~ 400 MeV, one has 
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In Fig. we show our results for fp calculated from 
substituting our solution for £(fc 2 ) into eq. I|3.8I) . In the 
walking limit, fp has been shown to satisfy a relation 
similar to eq. |3.6p . i.e., it is exponentially smaller than 
the scale A. We display, as the dotted curve, the fit from 
Ref. jg for the walking interval 0.89 < a* < 1.0, given by 
eq. (|3.6() with c = 1.5. Our results show the change from 
this walking type of behavior as a* increases above 1.2; as 
a* increases from 1.0 to 2.5, fp/A increases substantially, 
from about 3 x 10~ 3 to about 0.08. This is similar to 
the factor by which we found that EA increased as a» 
increased through this interval. 

The strong increase in E/A and fp/A. as a* ascends 
from the value 0.89 near the walking limit to the value 2.5 
deeper within the confinement phase is easily understood 
as reflecting the removal of the extreme exponential sup- 
pression evident in eq. 1)3. (jf) and its analogue for fp for 
a* — a cr — > + . One does not expect such a dramatic 
change in the ratio E//p over this interval, and this ex- 
pectation is borne out by our calculations. In Fig. [3] we 
show the ratio of E//p, which increases gradually from 
about 2.6 to 3.9. The fact that we find a ratio comparable 



FIG. 3: Plot of the ratio E//p for 1 < a* < 2.5. 



to the observed one in actual QCD, given by eq. (13. 711 . 
can be understood as a consequence of the property that 
much of the strong dependence on Nf divides out in this 
ratio. 



IV. CALCULATION OF MESON MASSES VIA 
THE BETHE-SALPETER EQUATION 

A. General Discussion 

We denote the wavefunction for a hadron with a given 
flavor combination for the generalized tt, p, etc. as 
follows. Define the flavor vector f a = (f a \ f aN '). 
Recall that in the confined phase the global symmetry 
SV(Nf) L x S\J(Nf)R x U(l)v is broken spontaneously 
to SU(Af/)y x U(l)y. We drop the explicit subscript V_ 
on SV(N f ) v henceforth. With regard to SU(N f ), a // 
meson with a given J p (where J denotes the spin, and 
P and C are the parity and charge conjugate eigenval- 
ues) is described via the Clebsch-Gordan decomposition 
Nt x Nf — 1 + Adj, where 1 and Adj denote the singlet 
and adjoint representations. 

Let the generators of the group SU(A/) have the 
standard normalization Tr(TjTj) = (1/2)5^-. Then the 
hadrons transforming according to the adjoint repre- 
sentation of SU(Af/) are comprised of (i) the set of 
Nf(N f - 1) states 



hv-ij — 



1 



/Nr. 



No 
a=l 



(4.1) 



where is the Nf x Nf matrix with a 1 in the i'th 
column and j'th row, with 1 < i, j • < Nf, i ^ j, and T 
specifies the type of particle (pseudoscalar, vector, axial- 
vector, scalar), and (ii) the Nf — 1 states corresponds to 
the generators of the Cartan subalgebra of SU(iVj) given 
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TABLE I: Data on relevant qq mesons whose masses are compared 
with Bethe-Salpeter calculations. n 2S+1 Lj is standard spectro- 
scopic notation, where n denotes radial quantum number. The 
symbols adj. and sing, denote the adjoint and singlet represen- 
tations of the SU(2) V- isospin flavor symmetry group. Masses are 
given in MeV, from |4.'J . The last column lists the mass divided by 
a typical hadronic scale, f n . 



n 2S+1 Lj 


J PC 


RSU(2) V 


name 


M 


M/U 


l 3 Si 


1" 


adj. 


P 


775.8 ±0.5 


8.40 


l 3 Si 


1"" 


sing. 




782.6 ±0.1 


8.47 


I 1 Pi 


1+- 


adj. 


bi 


1229.5 ± 3.2 


13.3 


I 1 Pi 


1+- 


sing. 


hi 


1170 ±20 


12.7 ±0.2 


l 3 Po 


0++ 


adj. 


do 


984.7 ± 1.2 


10.7 


l a Po 


0++ 


sing. 


fo 


~ 600_ 200 


6 5+ 6 - 5 

D -°-4.3 


I 3 Pi 


1++ 


adj. 


dl 


1230 ± 40 


13.3 ± 0.4 


I 3 Pi 


1++ 


sing. 


fl 


1281.8 ± 0.6 


13.9 



by the traceless Nf x Nf matrices 

T dk = [2k(k + l)]- 1 / 2 diag(l, 1, .., 1, -k, 0, ..,0) (4.2) 

where there are k l's and 1 < k < Nf — 
1. That is, T dl = (l/2)diag(l,-l,...,0), T d2 = 
(2v / 3)" 1 diag(l,l,-2,0,...,0), etc. Because of the 
SU(Nf) flavor symmetry, it does not matter which of 
these — 1 hadrons with a given T we use. We shall 
refer to these as the 7V/-generalized p, u>, etc. In partic- 
ular, the spectrum of mesons includes a set of Nf — 1 
Nambu-Goldstone bosons (NGB's) with L = S = 
and J PC = h , transforming according to the adjoint 
representation of SU(A^). The corresponding h sin- 
glet with respect to SU (Nf), i.e., the generalized rf , is 
not a Nambu-Goldstone boson because the correspond- 
ing U(\)a symmetry is anomalous. Our analysis of me- 
son masses is for the lowest-lying // states. In future 
work one could also consider radial excitations, Regge 
recurrences, pure gluonic states (glueballs) and the gen- 
eral coupled situation in which glueballs and // mesons 
of the same J PC mix. 

In QCD, there are several (light-quark) qq mesons that 
are of interest here. For the reader's convenience, we list 
these in Table [I] A notation for the various states in 
the case of general Nf massless quarks is Sr, Pr, Vr, 
and Ar, standing for "scalar, pseudoscalar, vector, and 
axial- vector" , where the subscript R denotes the repre- 
sentation - adjoint or singlet - under the SU(iV/) flavor 
symmetry group. The experimental and theoretical sit- 
uation concerning the ++ isoscalar meson /o has been 
the subject of much discussion over the years; indeed, 
this state may involve mixing with qqqq mesons [45^ . Be- 
cause of the complications in the analysis of this state, 
and the expected complications in a realistic analysis of 
its ^/-generalization, the SU(iV/-)-singlet ++ meson, we 
do not attempt to treat this in our current study. 

As will be seen below, in the Bethe-Salpeter equation 



that we use to calculate the masses of the mesons, the 
flavor-dependent structure is simply a prefactor. Hence, 
the solutions of this equation have the property that, for 
a given radial quantum number and spectroscopic form 
2S+1 Lj, the SU (Nf) flavor-singlet and flavor-adjoint 
mesons are degenerate: 

M(n 2S+1 Lj; flav. adjoint) = M(n 2S+1 Lj; flav. singlet) 

.(4.3) 

In view of this, we henceforth drop the subscript R 
and simply write V rather than Vfiav.adj. or Vfiav.sing., 
etc. Note that this is different from the prediction from 
SU (Nf) flavor symmetry (with degenerate quarks and 
electroweak interactions turned off) that the members of 
a given representation of SU (Nf) are degenerate. Exper- 
imentally, except for the pseudoscalar mesons, the light- 
quark isospin-adjoint and isospin-singlet qq mesons are 
nearly degenerate. The physical u> meson is very nearly 
a singlet under isospin SU(2), so a measure of this pre- 
dicted degeneracy for the ground state 1" mesons is 
(M UJ -M p )/[(l/2)((M UJ +M p )] = 0.87xl(T 2 , quite small. 
Similarly, (M fl -M ai )/[(l/2)((M fl +M 01 )] = 0.04±0.03 
and (M hl -M bl )/[(l/2)((M hl + M bl )} = -0.05±0.02. So 
for these states the prediction from our Bethe-Salpeter 
technique for the special case Nf = 2 massless quarks 
is in agreement with the data for light-quark mesons in 
QCD. 

The situation with the ^ mesons is quite different. 
Since the SU(Nf) flavor-singlet mesons are not NGB's, 
owing to the anomalous nature of the U(1)a symmetry, 
they are split by a large mass difference from the flavor- 
adjoint NGB's. In this noted above, the semi- 
perturbative Bethe-Salpeter analysis does not contain the 
relevant physics involving instantons, and hence its pre- 
diction is far from reality. For this reason we do not 
consider the flavor-singlet h mesons here. As regards 
the flavor-adjoint mesons, since we assume massless 
fermions and have turned off electroweak interactions, 
the mass Mp of the flavor-adjoint pseudoscalar mesons 
is exactly zero in our calculations. 

The pion decay constant f w provides a convenient mass 
scale with which to normalize the hadron masses. For 
comparison with our results calculated in the case of gen- 
eral larger Nf, we list in Table [I] the masses of the qq 
mesons divided by For later use we also record the 
ratio 



= 1.59 ±0.05 



(4.4) 



This is slightly larger than the prediction M ai /M p = 
V2 ~ 1.414 from a combination of vector meson domi- 
nance and spectral function sum rules ^(|. Also, 



M a 



= 1.27 . 



(4.5) 



An interesting result of the calculations of meson 
masses in the walking limit in Ref. |8j was that the ratios 
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of these masses to fp are rather constant. Specifically, it 
was found that in for 0.89 < a* < 1.0, 



and 



so that 



and 



My 
f P 

Ma 
fp ' 

Ms 
fp 

Ma 
M v 

Ms 
My " 
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11.5 , 



4.1 , 



1.04 



0.36 



(4.6) 



(4.7) 



(4.8) 



(4.9) 



(4.10) 



where the theoretical uncertainty is several per cent. 
These ratios may be compared with the values in regular 
QCD which, as far as the light-meson spectrum is con- 
cerned, are close to the values that they would have in the 
Nf = 2 chiral limit (with the understanding that the pion 
masses would actually vanish in this limit if electroweak 
interactions are turned off, as assumed here). For the 
purpose of this comparison, we do not try to use the in- 
ferred chiral-limit value of f^ ^(jj since to be consistent 
we would have to do the same for the mesons themselves. 
For the comparison between the extreme walking limit 
(WL) and QCD, we have 



(M v /fp) WL 
(M P /U) 



(Ma/ fp) 



WL 



and 



(M ai /U) 



(M S /fp)wL 

(M aa /U) 



1.3 



0.86 



0.38 . 



(4.13) 



A major output of the present work is the elucidation 
of how, as Nf decreases and a* increases, the ratios of 
meson masses to fp begin to shift toward their QCD 
values. 



V. CALCULATIONS 

Next, we describe the details of our solution of the 
Bethe-Salpeter equation and the resulting masses of // 
mesons. 



A. Bethe-Salpeter amplitudes 

We introduce the Bethe-Salpeter amplitudes \ f° r the 
scalar (S), pseudoscalar (P), vector (V), and axial- vector 
(^4) bound states of quark and anti-quark as follows : 



(Q\T1> afi (x+) ^' J (x^ \S a (q)) 
S'-mi-x f d4 P _-i 



(5.1) 



V2S{(T a y f e 



(2tt)< 



(0| T ip a fi(x+) V^Or-) \P a (q)) 

V2si(T a y; e-^- x 1 ''''>' 



(5.2) 



(2tt)< 



e w ' r [x(p)(p;g)] Q /3 



T^ afl (x+) jp(x.) \V a (q,e)) 
d 4 p 



(2tt)< 



(5.3) 

vp ' r \X(v)(p;q,e)] a /3 , 



T^ afi (x+) i;p{x-) \A a (q,e)) 

d P —ip-r 



- \/2S j (T ) f> p- iq X ( — 
-V2 dl {l a ) f e J w 



(5.4) 

[X(A)(p;<?,e)]a/3 , 



where x± = X ± r/2, and (a, (3), (/, /'), and (i, j) 
denote the spinor, flavor, and color indices, respectively. 
A a represents flavor structure of the bound states. In the 
case of flavor-adjoint bound states, T a is the generator of 
SU(Nf), while in the case of flavor singlet bound states, 

(T a ) f ! is the identity 1. 

We expand the BS amplitude x m terms of the bispinor 
bases T" 1 and the invariant amplitudes as follows : 



(4.ii) [x{s,p)ip;q)) a p = Y.\ r \s,p)(P^)] ap x { (lp)&<i) 



(5.5) 



(4- 12 ) [X(v,A)(p;q,e)] a0 = £ [i^ (p; <?, , 



a/3 



HV,A) 



(p; ?)• 

(5.6) 

The bispinor bases can be determined from the spin, 
parity, and charge conjugation properties of the bound 
states. The explicit forms of F\s)> ^\p)> ^(v)' an< ^ ^l A ) 
are summarized in Appendix lAl 

We take the rest frame of the bound state as a frame 
of reference: 



q» = (M B , 0,0,0), 



(5.7) 



where Mb represents the bound state mass, i.e., 
Ms, Mp, My and Ma for scalar, pseudoscalar, vec- 
tor and axial- vector meson masses, respectively. After 
a Wick rotation, we parametrize p^ by the real variables 
14 and x as 



p ■ q = iM B u , p 



2 2 
-U — X . 



(5.8) 
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m a,v,s,p/a 



P-Y2 



FIG. 4: A graphical representation of the HBS equation in 
the (improved) ladder approximation 



Consequently, the invariant amplitudes x^ can be ex- 
pressed as functions of the variables u and x: 



(0 _ v (») c \ 

X(s,p,y,A) — ^{s,p,v,a)\ u ^ x >- 



(5.9) 



From the charge conjugation properties for the BS am- 
plitude x an d the bispinor bases defined in Appendix lAl 
the invariant amplitudes x^ ( M : x ) are shown to satisfy 
the following relation: 



(i) i \ (i) / \ 

X(s,p,v,a)( u > x ) — X.(s,py,A)\ u -> x > 



(i) 



(5.10) 



B. Homogeneous Bethe-Salpeter equation 

The homogeneous Bethe-Salpeter (HBS) equation is 
the self-consistent equation for the Bethe-Salpeter am- 
plitude, and it is expressed as (see Fig. 0} 

T X = K X - (5.11) 

The kinetic part T is given by 

T(p;q) = Sj 1 ( P + q/2)(g>Sj\p-q/2) , (5.12) 

where the BS kernel K in the improved ladder approxi- 
mation is expressed as 



{p-k) 2 V ( P -k) 2 

In the above expressions we used the tensor product no- 
tation 



{A®B) X = A X B 
and the inner product notation 
d A k 



Kx (p; q) 



K(p,k) x(k;q) 



(5.14) 



(5.15) 



It should be noted that the fermion propagators included 
in T in eq. Ij5.12|) have complex-valued arguments after 
the Wick rotation [47j . The arguments of the mass func- 
tions appearing in the two legs of the Bethe-Salpeter am- 
plitude are expressed as 



(p±q/2) 2 



x '-[^r) TiuM B . (5.16) 



0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 



M A /A 

My/A •••• + ••• 

Ms/A — B- 

- Mp/A X 




FIG. 5: Values of meson masses divided by A calculated from 
the Schwinger-Dyson and Bethe-Salpeter equations. 



In general, it is difficult to obtain mass functions for com- 
plex arguments by solving the Schwinger-Dyson equation 
in the complex plane, especially because of the analytic 
structure of the running coupling in the complex momen- 
tum plane. However, in the case of large Nf QCD, the 
analyticity of the two-loop running coupling constant |48| 
makes it possible to solve for the mass function in the 
complex plane. This leads to the following approxima- 
tion, in accordance with eq. t\'d. 5|) 8J: 



Re [a(Xe l29 )] = a, 8(A 2 
Im[a(AV 2e )] = 0. 



X), 



(5.17) 
(5.18) 



where no confusion should result from the use of the sym- 
bol 8 on the right-hand side of eq. H5.17fl as the step 
function. . 



C. Numerical results 

We next present the results of the numerical calcula- 
tions for the masses of the mesons. We solve the homoge- 
neous Bethe-Salpeter equation as an eigenvalue problem, 
namely, the Bethe-Salpeter amplitude as an eigenfunc- 
tion and the mass of a bound state as an eigenvalue, 
denoted generically as Mb- Because the eigenvalue M B 
appears nonlinearly in the equation, we use so-called fic- 
titious eigenvalue method |23 to obtain the value of Mb ■ 
For details of numerical method to solve HBS equation, 
see Ref. 0] . In Fig.0 we show the values of meson masses 
divided by A calculated from the Schwinger-Dyson and 
Bethe-Salpeter equations in the range 0.9 < a* < 2.5. 
In Fig. H3 we plot the values of Mb / fp in the range of 
0.9 < a* < 2.5. In Fig. [7] we plot the meson mass ratios 
M A /M V and M s /M v in the range of 0.9 < a* < 2.5. 
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M A ,V,s/f P 



M A /f P o 

My/fp + 

Ms/fp □ 




-Q-H- 



1.2 1.4 1.6 



FIG. 6: Values of meson masses divided by Fp calculated 
from the Schwinger-Dyson and Bethe-Salpeter equations. 

M A ,s/M Y 

i i i i i i i i i 

1.4 ~ , 

M A /M V 

M S / My 



thereby approaching to within about 35 % of the 
value 10.7 in QCD for M a Jf^. 

3. My/fp decreases from about 11 to 9, rather close 
to the value 8.5 for M p j '/„■ and M^/fn in QCD. 
As is evident from Fig. H3 this ratio My/ f P is 
roughly constant in the upper end of the interval 
of a* values that we study. 

4. Ma/ fp behaves non-monotonically, first decreas- 
ing from roughly 11.5 to 10, but then increasing to 
about 11, within about 20 % of the average of the 
values in QCD for the isospin-triplet and isospin- 
singlet axial- vector mesons, 13 for M ai / ' f-w and 14 
for M fl /f n . 

5. Thus, the ratios M A /M V and M s /M v , which were 
found in Ref. to have the respective values 1.04 
and 0.36 in the walking limit, both increase in the 
interval of a* that we study, reaching about 1.17 
and 0.74, respectively, at a* = 2.5. For compari- 
son, these ratios are approximately 1.6 and 1.3 in 
QCD (cf. eqs. (|4.4() and (|4.5l0 . Although the value 
of the ratio Ms /My at a* = 2.5 is farther from 
its QCD value than is the case with M^/My, it 
is increasing somewhat more rapidly as a function 
of cv*, consistent with eventually approaching the 
QCD value. 



VI. CONCLUSIONS 



1.2 1.4 1.6 1.8 2 2.2 2.4 



a* 



FIG. 7: Ratios of meson masses calculated from the 
Schwinger-Dyson and Bethe-Salpeter equations. 



Our calculations yield a number of interesting results. 
We summarize these for the changes in these meson 
masses as a* increases from 0.9 to 2.5 as follows. 

1. The ratios of the meson masses divided by A in- 
crease dramatically, by factors of order 10 2 , ap- 
proaching values of order unity at a* = 2.5. This 
amounts to the removal of the exponential suppres- 
sion of these masses which had described the walk- 
ing limit at the boundary of the non-Abelian phase, 
as one moves away from this limit into the interior 
of the confined phase. 

2. Ms/fp increases monotonically from about 4 to 7, 



In this paper we have considered a vectorial S\J(N C ) 
gauge theory with Nf massless fermions transforming 
according to the fundamental representation and have 
studied the shift in behavior from walking that occurs in 
the region near the boundary between the confinement 
phase and the non-Abelian Coulomb phase to the QCD- 
like behavior with a non-walking coupling. Specifically, 
we have used the Schwinger-Dyson and Bethe-Salpeter 
equations to calculate the dynamically induced fermion 
mass E, the spontaneous chiral symmetry breaking pa- 
rameter fp, and the masses of the lowest- lying qq vec- 
tor, axial-vector, and flavor-adjoint scalar mesons. We 
have investigated how these change as one decreases Nf 
below Nf^ cr , or equivalently, increases a* above a cr , to 
move away from the above-mentioned boundary into the 
interior of the confinement phase. Our results show the 
crossover between walking and non-walking behavior in 
a gauge theory. 

There are a number of interesting topics for future re- 
search using the methods of this paper. It would be use- 
ful to construct a kernel for the Bethe-Salpeter equation 
that could include more of the relevant physics, including 
instantons effects. Work is underway on this. It would 
also be worthwhile to calculate the masses of radially 
excited mesons and mesons with internal orbital angu- 
lar momenta L > 2, as well as glueballs and the mixing 
between glueballs and qq mesons. We anticipate that 
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the results of these calculations would exhibit the same 
general properties that we have found with the ground- 
state qq mesons, but it would be interesting to confirm 
this expectation explicitly. Another project would be to 
connect our study of the region in Nf where there is a 
crossover from walking to nonwalking behavior, to the 
region around Nf — 2. However, when one moves this 
far away from the walking regime, one loses a simplifying 
features of our calculation, namely the fact that we do 
not have to use an infrared cutoff on a. Given that lat- 
tice gauge theory methods provide an ab initio framework 
for calculating hadron masses, we hope that the lattice 
community will extend early efforts such as those of Rcf. 

and carry out a definitive study of hadron masses 
in QCD with an arbitrary number of flavors. It would 
be of considerable interest to compare the results of the 
lattice calculations with those obtained from solutions of 
Schwinger-Dyson and Bethe-Salpeter equations. 
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APPENDIX A: BISPINOR BASES FOR SCALAR, 
PSEUDOSCALAR, VECTOR, AND 
AXIAL- VECTOR BOUND STATES 

In this appendix we show the explicit forms of the 
bispinor bases for the scalar, pseudoscalar, vector, and 



axial-vector bound states. Here we use the notation 
% = q^l^B with Mb being the mass of the bound 
states, and [a, 6, c] = a[6, c] + 6[c, a] + c[a, b]. 

Bispinor base for the scalar bound state ( J = ++ ) 
is given by 

T(S) = 1, r^=ji, Tf s) =i(p-q), 

rfs) = \ \Ul (Ai) 

and that for the pseudoscalar bound state (J PC = 0~ + ) 
is given by 

T(p) = 75, I^ P) =i> {p-q) 75, Tf P} = I 75, 
rfp) = i[M75. (A2) 

Furthermore, for the vector bound state (J PC = 1 ) 
we use 

rf F) = i^^j],rf y) = (e- P ), T\ v) =i>{e-p), 

r 7 iv) =i^-q)(e-p), Tf v) = ^J}(e-p), (A3) 
and for the axial-vector bound state (J PC = 1 ++ ) 

rU) = t n> r U = ^hs, rf A) = (p-q) 75, 

r (A) = 3j [/'^' I! 75, = {e-p) (p- q) 75, 

r (A) = i>( e ■ P) 75, t\ A ) = i (e -p) (p • 9) 75, 

rf A) = 5^|](e-p) (p-g)75- (A4) 
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